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Abstract 

The behavior of a spheroidal vesicle, in a plane shear flow bounded from 
one side by a wall, is analysed when the distance from the wall is much larger 
than the spheroid radius. It is found that tank treading motions produce a 
transverse drift away from the wall, proportional to the spheroid eccentricity 
and the inverse square of the distance from the wall. This drift is independent 
of inertia, and is completely determined by the characteristics of the vesicle 
membrane. The relative strength of the contribution to drift from tank-treading 
motions and from the presence of inertial corrections, is discussed. 
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I. Introduction 

An important topic in the study of suspensions is to understand the ability 
of shear to induce particle transport perpendicular to the velocity lines |2|. 
One of the mechanisms responsible for this kind of transport is the presence 
of lift forces due to inertia, pushing the particles away from the walls || (for 
recent references, see 0). It is well known that the symmetry properties of the 
Stokes equation, which governs the dynamics of a purely viscous fluid, do not 
allow for the presence of lift perpendicular to the flow lines, at least for the case 
of spherical particles suspended in a plane shear flow However, mechanisms 
different from inertia may allow for the symmetry breaking, necessary for the 
production of a transverse drift. 

A mechanism that is important, if one deals with a suspension of deformable 
objects, is the ability of shear to induce, under appropriate circumstances, a 
fixed orientation of the particles in suspension. If the deformable objects are 
vesicles filled with some other fluid, like in the case of blood, red cells, a fixed 
orientation can be attained when the vesicles are in a state of tank-treading 
motion ||. In this state, the membrane and the fluid inside circulate in a 
steady fashion around the cell interior, while the cell shape and orientation 
remain constant. 

Taking a comoving cartesian system {x\, X2, £3}, with origin at the cell cen- 
tre, and oriented in such a way that the unperturbed velocity field has compo- 
nents: v = (0, 0, aX2) : there will be a velocity perturbation v such as, if the cell 
is exactly spherical and inertia is neglible: ^1,2(^3) = —vi.2{—x^). Under these 
conditions, a plane wall perpendicular to the xi axis, will produce a correction 
v 1 to the perturbation field, such as, by symmetry v\ = at the cell centre. An 
ellipsoidal cell with fixed orientation, instead, will lead to symmetry breaking. 
Hence, there will be a drift f L = 1^(0), which will depend on the cell eccentricity, 
on the boundary conditions at the cell surface and, of course, on the distance 
from the wall I. Notice however that a non spherical shape is not by itself a 
sufficient condition for lift. A rigid non-spherical object for instance, will not 
keep an orientation fixed in time, carrying on instead a kind of flipping motion 
, whose contribution to lift is going to be much reduced. 
Inertia does not play any role in the mechanism for lift outlined above, and 
it will appear that its effect can be disregarded, when the particle is sufficiently 
close to the wall. Actually, depending on the circumstances, the contribution to 
drift from tank treading motions may even result to be dominant. The drift of 
a spherical particle in a bounded shear flow, is purely due to inertia; in the case 
of a particle in a Couette gap with thickness L: v^(l) = f(l/L)Re p aR, where R 
is the particle radius, Re p — 9L ^-, with v the kinematic viscosity of the external 
fluid, is the particle Reynolds number, and f(l/L) is at most of the order of a 
few tenths Q. In the case of a tank-treading motion generated drift, one has 
instead, in place of Re p , some parameter describing the non-sphericity of the 
particle, which, like in the case of red cells H, is not necessarily small. 
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In this paper, the simplest possible case of a neutrally buoyant, almost spher- 
ical vesicle in a purely Newtonian solvent, is taken into consideration. This 
allows a perturbative analysis around the standard case of a spherical particle. 
The vesicle is taken at a distance I from the wall such as: R <C I -C y/v/et', in 
this range, it is possible at the same time to neglect inertia, and to treat the 
effect of the wall, as a correction to the velocity perturbation due to the vesicle. 
The technique is similar to the one adopted by Ho and Leal in the case of 
a spherical particle in a bounded shear flow. 

The full problem of the determination of the spheroid shape under the com- 
bined effect of the stresses in the solvent, the membrane and the fluid inside 
is not treated (for reference about this problem, see e.g. Jll], |l2)). Rather, an 
axisymmetric, ellipsoidal vesicle shape is assumed. Thus, no information on the 
nature of the membrane and of the fluid inside is utilised. However, two kinds 
of motions at the vesicle boundary are analysed: one which is area preserving, 
mimicking the behavior of an inextensible membrane, like that of blood cells, 
and one in which points at the surface move with uniform angular velocity, 
simulating the case of an immaterial interface. 

II Stokes equation in vector spherical harmonics 

In a perturbative analysis around a spherically symmetric situation, it is 
worth while expanding the velocity field in terms of vector spherical harmonics: 

v(x) =v s (x)+v c (x)+v m (x) 
= Y,M m (x)xYim(e x ) + vf m (x)VY lm (e x ) + * V]Y* ro (e x )], (1) 

lm 

where Yi m (& x ) are standard, normalized spherical harmonics and e x = x/cc; the 
superscripts {sem} stand for scalar, electric and magnetic and come from the 
origin of this basis as a tool in the study of electromagnetic waves |l3| . 

This basis is clearly orthogonal and the components can be obtained in 
the standard way. Alternatively, one can write: 

< lm|x- v > 

v lm = ~2 > ( 2 ) 

vL = ~mTT) < Hv± ' (v ~ vS) > (3) 

and 

< lm\V ■ [x x v m ] > 

lm 1(1 + 1) ' () 

where the bra-ket notation < lm\f >= fi m (x) — f dn x Y l ^ n (e x ) f (x) , with dQ x 
the solid angle differential, is used, and is the angular part of the gradient. 
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At stationarity, an incompressible fluid in creeping flow conditions, obeys 
the time independent linearized version of the vorticity equation: 

V 2 [Vxv]=0 (5) 

together with the continuity equation: 

V-v = 0. (6) 

In terms of W; s ™ components, the vorticity equation reads: 



f^> = (7) 



f d 2 2 d Z(7 + 1) ^(1,2) 
V dx 2 x dx x 2 



where: 



f^=-vL + ^f m and f^=-vZ, (8) 
while the continuity equation takes the form: 



x^ + 3vf m - l l±^vt m = (9) 
da; x z 

From Eqns. (7-9), one obtains the "outside" and "inside" solutions: 

Q*lm x blm*E ? 

1(1+1) a lm% — J+i x > (10) 
Clm x ? 

and: 

J-2 



v lm 


= a 'lm X ' 




v lm 


1+3 
1(1+1) 




V lm 


= C lm X 





,1+2 _ b l mx l 



I 



(11) 



The expression for the velocity perturbation by a spherical particle in a strain 
flow, is thus obtained. Remembering the correspondence between fully symmet- 
ric, traceless Z-tensors and spherical harmonics with given Z, the contraction: 
X GL , t j X % X j 2 where otij = otji and an — (the convention of summation over re- 
peated indices is adopted) can be expressed as a linear combination of spherical 
harmonics with I = 2. One finds then from Eqn. (10): 

Vt = {^ ax 5 + T^bx^x^jkXjXk - ^bx~ 5 aijXj, (12) 

which is, for a — —1^2 the velocity perturbation due to a sphere of unitary 
radius put in a strain flow: Vi = ^aijXj [ fl4| . 
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III. Velocity perturbation by an ellipsoidal vesicle 

Boundary conditions at the wall 

An ellipsoidal cell in a shear flow will feel the effect of the strain component 
of the flow, which will tend to align strain and ellipsoid axes, while the vortical 
component will tend to make the cell rotate. If the viscosity of the fluid inside 
or the rigidity of the membrane are too large, the vesicle will tend to rotate as 
a rigid body, otherwise, it will be in a state of tank-treading motion. 




V 

Figure 1: Cell orientation in a plane shear flow v = [0, 0, aX2\ 
In this configuration the symmetry axis of the ellipsoid will lie in the plane X2X3, 
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at an angle 9 with respect to the X2, with 9 — > ±7r/4 for vanishing resistance 
of the vesicle to deformation, and 9 — > ±7r/2 at the threshold for the transition 
to flipping motion In the two limits, the long side of the ellipsoid tends 
to be parallel, respectively, to the expanding direction of the strain and to the 
unperturbed velocity field. 

Taking a reference system {x'ix' 2 x' 3 } at the cell centre, as shown in Fig. 1, 
with x' 3 along the symmetry axis of the ellipsoid and x[ = X\, the equation for 
the cell surface can be written in the form: 

{l-e){x l2 +x' 2 ) + {l + e)xf = R 2 (13) 

For small values of the eccentricity e, the distance, from the origin of a point on 
the surface, with given elevation x' 3 (and corresponding angular coordinates in 
the unprimcd frame 9 and <f>), can be written in the form: 



x ( ^ fl (i + i(i-2(4) ; 



= R(l + i(l-2 (C03Sxs - R r Sx ^ )), (14) 

In a state of tank-treading motion, the membrane is assumed to move following 
two possible laws. In the first case, for each value of x±, points on the surface 
are taken to move at constant speed: 

v B = ^BaR(x 1 ), with: R{x x ) 2 = R 2 - (1 - e)x\, (15) 

remaining on the elliptic trajectory, which corresponds to the given section at 
constant x\ of the cell surface. Such motion is locally area preserving, and ap- 
proximates the behavior of an inextensible membrane, with negligible resistance 
to bending and shear. For small e, one expects B ~ 1 4- (3e; thus, for e = 0, the 
expression for vb reduces to that for a sphere immersed in the given shear flow. 
In terms of components: 

vb = vb(0, , X ' 3 , , X ' 2 ) (16) 



where 7 = \/ jxf ■ F° r small e 



v B * f [o, + W + 7£T + ^))»4(i + <P + i^T- ^3-))] ( 17 ) 



-± 



where x\ = x' 2 2 + x' 3 2 . Using Eqns. (2) and (17), one finds: 



— eax — — 



(18) 



G 



and: 

V B^ £ «(J-^)^3 (19) 

It will appear that only the / = 2 part of vb must be taken in consideration to 
calculate the lift. Using Eqn. (3), a simple calculation leads to the result: 

v B = v£ - jjeax^ - g[0, 4, x' 2 ] + H.H., (20) 

where H.H. indicates higher harmonics. In terms of unprimed variables and 
components: 

m 5 (x\ - 4) sin 20 + 2a:2^3cos20 
vb = v B - -ea X - 

cot - - - 

- — [0, -£ 2 sin26> + a;3COs26»,X3sin26> + 2;2COs26»] + H.H. (21) 

The second membrane behavior that is considered, is characterized by a constant 
angular velocity motion: 

Bet ex etc 

vb = — [o, -x s /7, 7^2] ~ - (1 + /?e) [o, -4, 4] - — [o, 4 , 4] , (22) 

In terms of unprimed variables: 

v B ~ |(l + /3e)[0,-x 3 ,x 2 ] 

cue - - - - 

— — [0, -x 2 sin 26 + x 3 cos 29, x 3 sin 29 + x 2 cos 29} . (23) 

The behavior described in Eqns. (22) and (23) should mimic the motion of a 
droplet; however, this is also the approximation for an area preserving mem- 
brane motion, used by Keller and Skalak in [||. Thus, comparing the velocity 
perturbation and the lift produced by the two membrane motions, prescribed 
in Eqns. (17) and (23), gives, among the other things, an idea on the accept- 
ability of such an approximation, also in view of the simplifications that become 
possible in the large e regime. 

Calculation of the perturbation in an infinite fluid 

If the membrane velocity -vb{9,4>) and the cell orientation 9, produced by 
the external flow v(#, <f>, x), are supposed known, the velocity perturbation v 
can be obtained from the boundary condition at the vesicle surface x — x(9, </>): 

w B {9, 0) = 9(9, <j>, x(9, 4>)) + v(0, <P, x{9, </>)). (24) 

For small values of e, it is possible to solve Eqn. (24) perturbatively; to 0(e): 

v™ (0, <f>) = v(6», 0, R) + v (0) (9, <f>, R) (25a) 
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vg> (9, 4>) = v(0, <j>, x{0, 0)) + v (t » (0, 0, 3(0, 0)) + v (1) (0, 0, i?), (25b) 

where, for x{6,<j>), the 0(e) correct expression, provided by Eqn. (14) can be 
used. In Eqn. (25a) v^', is the velocity at the surface of the sphere, with 
x (0) (6, 4>) = R, immersed in the same flow: 

vg>(M) = f[0,-x<°\x<°>]; (26) 

In this way, v <0) is the velocity perturbation by a sphere in a shear flow, given 
by Eqn. (12): 

(0) 5/R 4 R 2 \ax 2 x 3 x aR 4 
V = 2 t - X*) ~ ^ [ °' ^ X2] (2?) 

Expanding Eqn. (25b) in vector spherical harmonics and using Eqn. (10), allows 
then to calculate v (1) for x > x(0, <p). 

For the first membrane motion, described by Eqn. (21), one has: 

8x s 



where: 



V * = f ( (/? + 2^ )£ + W [ °' ~* 3 ' X2] + VC m 

5 ixl - in ) sin 29 + 2x 2 x 3 cos 29 

vc = eax — - 

12 R 2 

eo; - - - - 

- — [0,-32 sin 26» + x 3 cos 26>,x 3 sin 26» + x 2 cos 26>] + H.H. (29) 

In Eqns. (28) and (29): 

fa = g!( 1 _ 2 ( c0s ^3- 2 sinfa2)^ (2Q) 

and the identification x (0) = x is allowed at the order in e considered. It 
appears (details are given in the appendix) that the condition on the absence 
of external torque on the vesicle, corresponds to the I = 1 components of v m , 
being identically zero. This is achieved choosing appropriately the coefficient f3, 
and consequently, the value of the tank treading velocity. The present analysis 
is directed to calculate the lift at large distance from the wall. Hence, from Eqn. 
(10), the leading contribution is from vf =2 terms. Substituting Eqns. (14), (27) 
and (28) into Eqn. (25b), keeping only 0(e) terms, and neglecting v m and / > 2 
components, one obtains for v (1) : 



5aSx /X2X3X 1 w a 
< x = R = vc + —jj-{—j^- -[0,x 3 ,x 2 ]J. (30) 

Using Eqn. (29) together with Eqns. (2-3), one finds: 



v ( x l s R = [(z§ - a£) sin 29 + 2x 2 x 3 cos 29] (31) 
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and: 

(i )0 _ mR 2 v (j4 ~ 4) sin 29 + (2 cos 29 + 30)0:20:3 

x=R — x 2 

__( coste3 _ sinfe2 )2^^_|_ _ _[ 0)2 . 3j X2 ]j. (3 2 ) 

All terms in Eqns. (31) and (32) are clearly 2 = 2 with the exception of the 
last line of Eqn. (32), which is a combination of v e and v m terms. After some 
tedious calculations, expedited by the use of Eqn. (3), the following expression 
for the / = 2 part of vjj. 1 ^, is obtained: 

„(Dc p2^/ giDjg /5xj x\ 17a;|\ _/ 25 cos20\ 1-2X3 \ 

The leading term at large distance is, in the notation of Eqn. (10), the contri- 
bution a,2 m x~ 3 to v s . From Eqns. (31) and (33), using Eqn. (10), one obtains 
therefore the large distance result: 



v 



eai? 3 x / sin 29 t 5x'{ 2hx\ 45x§\ /25 5 cos 20 \ x 2 x 3 \ 
x^\~^~\~14~ + ~56 56/ V28 + 4 J^ 2 "/' ( ' 



The same identical calculations can be carried on for the second kind of mem- 
brane motion described by Eqn. (23), with the final result: 

m eai? 3 x / sin 29 ( 5x? 152:? 10a;l\ /25 5 cos 20 \ a; 2 ^3 \ ,„ N 
v(, "^(^(lT + ^4 1 -V)-(28 + ^)^)- (35) 

One question that comes natural at this point is whether one could extend the 
analysis to a non perturbative regime, expanding directly Eqn. (24) into vector 
spherical harmonics and then, after imposing a cut-off on I and m, solving 
numerically the resulting linear system in the unknowns ai m , bi m and c; m . This 
could be useful to analyse the behavior of vesicles of a general, strongly non 
spherical shape, but it is not the most efficient way of proceeding, also because 
the matrix associated with the system is strongly ill conditioned. Of course, 
in the case of an ellipsoid immersed in a shear flow, undergoing a linear tank- 
treading motion like the one of Eqn. (22), the theory of Keller and Skalak § 
can be utilised to obtain analytic expressions for the velocity field. In the more 
general case, it turns out that it is possible to modify the initial vector spherical 
harmonics basis, in such a way that the velocity field is obtained directly, without 
having to invert any ill conditioned matrices [l5| . 

IV. Calculation of lift 

The calculations in the previous section referred to an unbounded flow sit- 
uation. A wall at a large distance from the vesicle will cause a correction v 1 to 
the velocity perturbation v obtained in the previous section. This can be cal- 
culated using the boundary condition provided by making the correction at the 
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wall, equal to minus the velocity given by Eqns. (34) or (35). The higher order 
corrections are obtained then, using the value, alternately on the vesicle surface 
and on the wall, of the previous calculated correction as boundary condition, 
in a series of images and counter-images, analogous to those of electrostatics. 
If the distance from the wall is large, it is possible however to stop at the first 
image. 

The technique at this point is standard |l^, ^6|; one introduces scalar and 
vector potential <f> and A, such as: 

v T =V^ + VxA (36) 

where: 

V 2 4> = and V • A = 0. (37) 

The first of Eqn. (37) is a consequence of incompressibility, while the second 
is a gauge condition. From here, the vorticity equation takes the form of a 
biquadratic: 

V 2 V 2 A = (38) 

Taking the wall at x 2 — I positive, parallel to the Xi i3 axes, it is useful to Fourier 
transform all quantities with respect to X\ and x 3 : /(x) = J d ^ 2 ^ 2 3 e l ^ klXl+k3X3 ' ) 

x/(fci, £3,2:2). The gauge condition on A takes then the form: 

i 3 = -^ii + M, (39) 
k 3 k 3 

where the prime indicates derivative with respect to x 2 . The gauge is definitely 
fixed by requiring that A does not contain any potential contribution a such as 
V 2 a = 0. With this condition, solution of Eqn. (38) gives: 

A = A(fci, k 3 ){x 2 - I) exp(fc(a; 2 - I)), (40) 



where k = \Jk\ + fc| . The first of Eqn. (37), instead, gives for <j>: 

4> = $(ki,k 3 )exp(k(x 2 -l)). (41) 

Using Eqns. (36) and (39), the expression for the velocity correction becomes, 
in terms of Fourier components: 

v\ = -§M; + jLA% - ik 3 A 2 + iktf, 
v\ = + ^A' 2 + 4>>, (42) 

«3 = ikiA 2 - A\ + ik 3 <p, 

and imposing the boundary condition v^fci, k 3 , 1) = — v(fci, k 3 , 1), one finds: 

vi = -^A 2 - k$, (43) 



v 3 = Ai - ik 



3 <: 
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where v(fci, £3) = v(fci, £3, I). Solution of this system gives: 

i[—k 1 k 3 v 1 +2ik 3 kv2~klv 3 ] 



_ kik 3 v 1 —2ik 3 kv 2 ~(k +k 3 )v 3 (44) 



2fc 2 

a _ i(k3Vi—kiV 3 ) 
A 2 - 2k 



At this point, one hnds v 2 , from substitution of Eqn. (44) into (42); for x 2 = 0: 
u 2 (fci, &3, 0) = —[ikilvi + (1 + kl)v 2 + ik 3 lv 3 ] exp(-kl) (45) 
The lift is found from inverse Fourier transform at x\ = 2; 3 = 0: 1)^(1) 

d 2 fc 

(27T) 2 



J (2^)2 ^2(^1 j ^3, 0). Of all the contributions to v entering Eqn. (45), only those 



proportional to sin 29 in v (1) give a nonzero result. Writing in an explicit way: 



27T rl-K 



v L (l) = — — / fcdfc / xdi / dyj / d"0exp(— ikx cosip — kl) 
(2txY J J q J J 

x [klvi cos <p + (1 + fcZ)w 2 + ifciu3 sin y>] (46) 



where i = \J x\ + x\ : x\ — xcosd, k\ — kcos(p and ip = 9 — (p. In Eqn. (46), 
one can write, from Eqns. (34) and (35): 

v= eaR3 f n 20 (2C1 2 + 2(D-C)x 2 cos 26 - Dx) [x cos 6, 1, x sin 9] , (47) 

x b 

where all terms giving zero contribution in the integrals of Eqn. (46) are disre- 
garded. The integrals in Eqn. (46) can be carried on analytically exploiting the 
following properties of the Bessel functions J„ (x) |l7j : 

/■27T 

/ dip exp(— i a cosip) — 27rJ (a) (48) 
Jo 

p2rr 

/ dip cosip &Kp(—ia cosip) — 27riJi(a); (49) 
Jo 

and: 



(J a 2 + 2 - a) v 
x»dxJ v (Px)exp(-ax) = (-1)*— ^ ^ J . (50) 

The final result is: 

eaCR 3 sin 2(9 



(51) 

Thus, the magnitude and the direction of the lift depend on the coefficient of 
the x\ terms in the expression for v (1) given by Eqns. (34) and (35). It is known 
||, and a simple argument is given in the appendix, that for sufficiently small 
values of e, 9 — ±~, with the plus sign when the ellipsoid is oblate, (e > 0) and 
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minus when this is prolate (e < 0). This leads then, both in the inextensible 
membrane case of Eqn. (34), and in the "droplet" case of Eqn. (35), a drift 
away from the wall, which is directly proportional to the inverse square of the 
distance from the wall, and to the eccentricity e: 

3£* (52) 
The values of the constant C are obtained form Eqns. (34-35); one finds in the 

25 j n — 15 
112 anCL ° - 28 ' 

V. Conclusions 



two cases, respectively: C = and C — — 



The main result that has been obtained in this paper, is that tank-treading 
motions are able to produce a transverse drift of vesicles in a sheared suspension. 
This effect is mainly localized near the walls, where it can dominate that of 
the inertial corrections, on which current theories on the lift of particles in 
suspension are based. It is possible to estimate the thickness of the region 
where this happens. Using the expression for the inertial drift of a spherical 
particle, in a bounded shear flow with < I < L: f(l/L) ~ f° (l/L — 0.5), 
where f° is a constant of the order of a few tenth M, one finds from Eqn. (51), 
that the drift from tank treading motions remains dominant as long as: 

I / CeL \i .„ . 

<(^-s) • (53) 



R \f°Re p R 

Thus, thanks to the smallness in most situations, of the particle Reynolds num- 
ber Re p , the width of the region where the drift from tank-treading motions is 
dominant can be so large that both approximations of large l/R and small e, 
used in this paper, can be satisfied at the same time. 

However, the most interesting situation is that of strongly non spherical vesi- 
cles, for which e = 0(1). It is known that, for large values of e, the vesicle tends 
to have its major axis aligned with the unperturbed flow (which corresponds to 
sin 29 = 0), and then to make the transition to flipping motion ||. Thus, there 
must be some critical e, for which, fixed the other parameters, the drift velocity 
achieves its maximum value. 

In the present theory, due to the smallness of e, one has: 6 ~ ±7r/4, and 
the effect of destruction of tank treading motions, produced by e becoming too 
large, is not accounted for. At the same time, the perturbative calculation 
leading to the large distance behavior of the velocity disturbance v, becomes 
unreliable. However, extrapolating Eqn. (53) to large values of e, the effect 
of tank-treading motions (when present) is expected to become dominant in 
all situations of suspensions flowing in narrow gaps. Thus a nonperturbative 
extension of the theory at large e would be advisable. 

Red cells are a physical system in which such a non perturbative theory 



could be applied. However, experimental observations |1S, 19 as well as theo 



retical analysis |p[ indicate that, due to their membrane viscosity and that of 
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the hemoglobin inside, these cells undergo tank treading motions only when im- 
mersed in very viscous solvents, or in the presence of very strong shear stresses. 
It is clear on the other hand, that inertia cannot account by itself for phenom- 
ena like the concentration of red cells near the axis of small blood vessels: the 
so called Fahraeus-Lindqvist effect Q . Using typical parameters for a red cell 
and a small vessel ||: R ~ 4/x, L ~ 100/i, a ~ 200sec _1 (i.e. v ~ lem/sec), one 
finds that a transverse drift of the order of L, would occur only after the cell has 
travelled already for several centimeters, much more than the length of a typi- 
cal small vessel. Thus, a fixed orientation and strong departure from sphericity 
remain the necessary ingredients for a sufficiently large lift, and some different 
mechanism for maintaining a fixed orientation in the cell must be found. 

The most important limitation of an analysis like the one in this paper, 
is the fact that the membrane motion and the cell shape are imposed from 
the outside. This, among the other things, makes impossible to analyse the 
behavior of the cell, when the distance from the wall becomes comparable to 
the radius R. At the same time, there remain open several questions on which 
kind of membrane motions and cell shapes are possible, for instance whether 
there could be membranes, whose response to an external shear field leads to 
a drift towards the nearest wall, instead of away from it. The two behaviors 
considered in this work lead to lifts of different magnitude, but both directed 
away from the wall. 

The fact that different membrane behaviors lead to the cell moving in one di- 
rection or another suggests the possibility of controlling the cell motion through 
the membrane stiffness and the internal fluid viscosity. This could be of great 
importance in microsurgery, where new methods to deliver drugs using micro- 
capsules that break and drop their content when they are on the selected target, 
are taken into consideration. Already now, it seems realistic to have cells whose 
stiffness varies with temperature and whose tendency to stay away from the 
walls could be controlled from the outside. Carefully designing the cell structure 
(and its response to external stresses or other stimuli), could allow in principle, 
a fuller control on the cell trajectory, without having to rely on internal motors 
and sources of energy. 
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Appendix: Torque and energy balance 

In a steady state of tank trading motion, a vesicle in a shear flow will feel no 
external torque, while the external work will be equal to the energy dissipated in 
the membrane and the cell interior. In this paper, the membrane is supposed to 
oppose no resistance to bending and to shear, and consequently, to be dissipation 
free. The simultaneous satisfaction of these two conditions fixes the magnitude 
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of the tank-treading velocity vb and the orientation angle 9. 

The external torque is calculated from the value of the stress tensor at 
a spherical surface enclosing the vesicle: 

Tij = fiiy(diVj + djVi), (Al) 

where fj, is the fluid density. The force and the torque on an infinitesimal el- 
ement dSi of the cell surface, oriented towards the outside, will be equal to, 
respectively: 

dll, = TijdSj and M t = ei jk XjdU k . (A2) 

From here, the total torque acting of the spherical surface in exam, can be 
expressed, after a few manipulations, in the form: 



(A3) 



M = [ivx J dtt x {{x x V](x • v) + [x x (x • V)v] - [x x v]}. 

The decomposition in vector spherical harmonics is used again: 

v = u s xY + v c VY + u m [x x V]F, (A4) 

where the notation: v scnl Y = ^2 lm v^Yi m is used. 

Substituting back into Eqn. (A3), one obtains the result that only the v m 
components contribute to M: 

M = - f ivx 5 (v m /xy J dO x V7 (A5) 

and, remembering the correspondence between spherical harmonics and irre- 
ducible tensors, it appears that only the I = 1 components survive. The no 
torque condition takes then the form, from Eqn. (10): 

cj=i = 0. (A6) 

Turning to the requirement of energy conservation and steady state, the equation 
for the balance between dissipation and external work can be written in the form: 



W 



= Jw B • d(n' - n) = o, (A7) 



where W is the membrane dissipation, taken equal to zero, and dll' is the force 
with which the fluid inside the cell acts on the infinitesimal element dS of the 
membrane. If the vesicle is perfectly spherical, it will rotate as a whole, so that 
W — 0, dll' = and Eqn. (A7) will coincide with the no torque condition 
M = 0. In fact, in this case only one parameter, vb, remains to be calculated. 

For small non-zero e, W is a linear functional of the external velocity field 
and can be decomposed into a contribution due to the strain S and another to 
the vorticity wofv: 

W = W s (9,e,S) + W ul (9,e,Lu). (A8) 
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If uj = the cell will tend to align with the strain axes of v, i.e.: 9 — 
therefore: 



W s (±^,e,S) = W s (9,0,S) = Q (A9) 

and the first nonzero derivative of W s is: d§d e W s (±j, 0, S). Turning to the 
vortical component, writing: 

J dSA M and W" = J dSA w , (AlO) 



M = / <LS i w and \V 

one observes that: 



A w = (c + 0(e))A M + 0(e 2 ). (All) 

Thus, if the no torque condition M = is satisfied, one has: W u (9, e, oj) = 
^d e 2 W(0,O,w) + C(e 3 ), and one finds from Eqn. (A7): 

Hence, almost spherical vesicles tend to carry on a tank-treading motion, with 
the symmetry axis alligned with one of the strain axes of the external flow. 
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